Astronomy <fc Astrophysics manuscript no. H4495 


February 2, 2008 


(DOI: will be inserted by hand later) 





Optical spectra of selected Chamaeleon I young 

stellar objects* 

fT) . Carlos Saffe 1 '**, Mercedes Gomez 1 '***, Sofia Randich 2 , Diego Mardones 3 , Paola 
O 



O 

(N 

m 



Caselli 2 , Paolo Persi 4 , and German Racca 1 

1 Observatorio Astronomico de Cordoba, Laprida 854, 5000 Cordoba, Argentina 
email: safTe@oac.uncor.edu, mercedes@oac.uncor.edu, german@oac.uncor.edu 

2 INAF-Osservatorio Astrofisico di Arcetri, L.go E. Fermi, 5 50125 Firenze, Italia 

email: randich@arcetri.astro.it, caselli@arcetri.astro.it 

> ■ 

\& ' 3 Departamento de Astronomi'a, Universidad de Chile, Casilla 36-D, Santiago, 

in 
o 
o. 

(^) , Istituto Astrofisica Spaziale e Fisica Cosmica, CNR, Via del Fosso del Cavaliere 

o 

, F00133, Roma, Italia 

email: persi@rm.iasf.cnr.it 

6 

H 



CO 

c3 



Chile 

email: mardones@das.uchile.cl 



Received Month XX, 2003; accepted Month XX, 2003 



Abstract. We present optical spectra of eight candidate brown dwarfs and a pre- 
viously known T Tauri star (Sz 33) of the Chamaeleon I dark cloud. We derived 
spectral types based on the strength of the TiO or VO absorption bands present 



in the spectra of these objects as well as on the PC3 index of 



Martin et al. 



1999). Photometric data from the liter ature are used to estimate the b olometric 



luminosities for these sources. We apply 



D'Antona fc Mazzitelli 



1997) pre-main 



sequence evolutionary tracks and isochrones to derive masses and ages. Based on 
the presence of Ha in emission, we confirm that most of the candidates are young 
objects. Our sample however includes two sources for which we can only provide 
upper limits for the emission in Ha; whereas these two objects are most likely 
foreground/background stars, higher resolution spectra are required to confirm 
their true nature. Among the likely cloud members, we detect one new sub-stellar 
object and three transition stellar/sub-stellar sources. 
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1. Introduction 



The field of star and planet formation has attracted much observational and theoretical 
interest, particularly during the last twenty years. Impressive technological advances 
have reached high spatial resolution and sensitivity limits and allowed us to observe true 
stellar embryos and their circumstellar environments (disks), where planets are supposed 
to form. Brown dwarfs (BDs) represent an intermediate class of objects between stars 
and planets. They have masses below the H- burning limit and then are unable to sustain 
stable nuclear fusion of H like stars. However, unlike planets, BDs once formed go through 
a slow gravitational contraction process that releases potential energy and makes them 
shine. Thus, their brightne ss decreases with time as the energy reservoir diminishes and 



Oppenheirmer et al. 



we adopt 80 Mj (1 Mj = 1(T 



they cool down. Following 
Mq) as the upper mass limit for BDs and 13 Mj as the lower mass limit. Solar metallicity 
objects wi th masses within the se limits can burn deuterium at least during part of the 



evolution (Bu rrows et al. 



1997T1 . Below 13 Mj there is no nuclear fusion of any element. 



Several invest i gations suggest a star -like formation mode for BDs (see, for example, 



Elmegreen 



1999; 



Natta fc Testi 



lar embryos lIReipurth fc Clarke 



I Pic kett et al. 



2001). However other two scenarios, ejection of stel- 
20(Jl| and planet-like formation in circumstellar disks 



2000) have been proposed. The current observational evi dence seems to 



provi de support to the stellar-like formation process (see, for example, 



Muench et al. 



20021) . although neither the embryos ejection nor the planet-like formation hypothesis 



can presently be disregarded on observational grounds. Understanding the formation 
process of BDs is very important as it can shed light on the mechanism of formation of 
both stars and planets. 

BDs are elusive because they are very faint. These objects are easier to detect when 



dark clouds (see, for example, 



Briceho et al. 


1998: 


Comeron et al. 


2000 


Briceho et al. 



2002). A complete census of these objects in a given star- forming region would allow 
us to carry out a meaningful estimate of their real number and thus a measurement of 
the IMF of the cloud. However a considerable observational effort is required to fully 
accomplish this aim. 

In this contribution we present results of optical spectroscopic observations of eight 
candidate young BDs in the Chamaeleon I dark cloud proposed by several recent infrared 
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surveys aimed at finding sub-stellar objects in this star-forming region (see, for example, 



Persi et al. 



2000; 



Gomez fc Kenvon 



2001). An optical spectrum of Sz 33, a previously 



known T Tauri star of the cloud, is also presented. Our final goal is to carry out a precise 
determination of the complete IMF of the cloud. 

In §2 we describe the observations and data reduction. In §3 we present our analysis 
a nd results. We d erive spectral types from the TiO and VO indexes and the PC3 index 
1 1999) which combined with published photometry allows us to place 



of 



Martin et al. 



the stars in the HR diagram. We conclude with a brief summary in §4. 



2. Observations and data reduction 

The observations were carried out from May to July 2001 in service mode with FORS1 
on the ESO-VLT UT1. This instrument has different configurations 1 . We used the MOS 
(multi-object spectroscopy) sub-mode in combination with the grisms GRISM— 600R 
and GRISM— 6001 to cover the spectral range between 5300 and 9000 A, approximately, 
on a 2048x2048 pixels CCD detector. These grisms provide resolving powers K— 1230 
and 1530, respectively. The corresponding spectral dispersions are 1.08 and 1.06 A/pix. 
Two order separation filters, GG435+31 and OG590+72, were employed with grisms 
GRISM— 600R and GRISM— 6001, respectively, to isolate the useful spectral range for 
each grism. 

MOS has 19 movable slitlets that can be displaced by linear guides to any position 
along the dispersion direction in the field of view. In addition, the instrument can be 
rotated around its optical axis allowing a wide variety of target configurations. The 
width of each individual slitlet is also adjustable. For our observations we selected a 
width of 1" for all the 19 slitlets. Even-numbered slitlets are 20" long and odd-numbered 
are 22" long. 

In addition to the VLT observations we obtained spectra of two relatively bright 
objects (ISO-Chal 126 and ISO-Chal 237, see Table 2) with the 2.15-m telescope at the 
CASLEO (San Juan, Argentina) in April 2000. We used the REOSC spectrograph in 
SD (simple dispersion) mode with a 300 1/mm grating blazed at 5000 A and a long 
slit opened to ~ 2". We registered the spectral range between 4500 and 8500 A on 
a 1024x1024 thinned TEK CCD detector at a dispersion of ~ 4 A/pix, resulting in a 
resolving power R= 750. 

We reduced the data using standard IRAF procedures 2 . We trimmed the CCD frames 
at each end of the slit, corrected for the bias level and divided by appropriate normalized 

1 See the FORS1+2 User Manual for a full description on the different observing modes of 

this instrument, available at: |http: / / www.eso.org/instruments / forsT/l 

2 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by 

the Association of Universities for Research in Astronomy, Inc. under contract to the National 
Science Foundation. 
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Table 1. Observed Objects 



Name 


a(2000.0) 


5(2000.0) 






IT R 


IT I 


Tel+Instru 


Other ID 




[h:m:sec] 


["'■■"] 






[sec] 


[sec] 






[GK2001] 8 


11 05 48.8 


- 76 40 17 






3050 


2950 


VLT+MOS 




ISO-Chal 126 


11 08 04.2 


- 77 38 43 






900 




CASLEO+REOSC 


[CCE98] 32 


[KG2001] 102 


11 09 49.3 


- 77 31 20 


2 


x 


3105 


2960 


VLT+MOS 




[SGR2003] 1 


11 09 53.1 


- 77 30 58 


2 


x 


3105 


2960 


VLT+MOS 




[GK2001] 30 


11 09 53.4 


- 77 28 36 


2 


X 


3105 


2960 


VLT+MOS 


ISO-Chal 220 


Sz 33 


11 09 53.9 


- 76 29 25 






3040 


2960 


VLT+MOS 


ISO-Chal 224, [KG2001] 110 


[GK2001] 31 


11 09 55.0 


- 76 31 12 






3040 


2960 


VLT+MOS 


ISO-Chal 225 


ISO-Chal 237 


11 10 11.9 


- 77 35 31 






900 




CASLEO+REOSC 


[CCE98] 48JOTS99] 45, [KG2001] 121 


[GK2001] 53 


11 14 23.6 


- 77 56 12 






3300 


3300 


VLT+MOS 





Note: We adopt the SIMBAD source denomination scheme which is based on the recommenda- 
tions of the IAU Commission 5 Task Group on Designations. 

flat-fields. The spectra were extracted using the NO AO task apextract with an aperture 
of 5 pixel radius. A sky subtraction was done by fitting a polynomial to the regions 
on either side of the aperture. Residual night sky lines were removed by interpolating 
across them. A non linear low order fit to the lines in the HeNeAr lamp (VLT spectra) 
or CuNeAr lamp (CASLEO data) was used to wavelength calibrate the spectra. Typical 
RMS for the wavelength solutions are 0.07 and 0.04 for the grisms GRISM— 600R and 
GRISM— 6001, respectively. We obtained values of ~ 0.5 for the RMS corresponding to 
the CASLEO solution. The VLT-FORS spectra have been flux-calibrated and normalized 
to the flux at 6800 A. CASLEO data are presented in count units as no flux standards 
were observed in this case. 



3. Data analysis and results 

3.1. The sample 



In Table 1 we list our sample stars together with their coordinates and observing logs; 
in Table 2 photometric data from the literature are provided. Seven of the nine stars in 
our sample were selected among young low m ass and brown dwarf candidates proposed 



by several recent infrared surveys of the clouc 


1 Cambresv et al. Ill998 


Persi et al. 2000 


Gomez & Kcnvon 


2001; 


Kenvon & Gomez 


2001 


). These objects, unlike most of the 



proposed candidates, have optical counterparts on the the Digitized Sky Survey 3 (DSS) 



3 Based on photographic data obtained using The UK Schmidt Telescope. The UK Schmidt 
Telescope was operated by the Royal Observatory Edinburgh, with funding from the UK Science 
and Engineering Research Council, until 1988 June, and thereafter by the Anglo- Australian 
Observatory. Original plate material is copyright (c) the Royal Observatory Edinburgh and 
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Table 2. Compiled Photometric Data and Hev Equivalent Width Measurements 



Name 


L a 


K b 


H b 


J b 


I c 


Ha (EW[A]) 


Na I 8183-8195 


[GK2001] 8 




13.94 


14.67 


15.22 




< - 5 


4.9 


ISO-Chal 126 




8.30 


9.67 


11.51 


14.39 


- 97.1 




[KG2001] 102 


10.62 


11.07 


11.45 


12.47 




- 87.2 


2.7 


[SGR2003] 1 




14.87 


15.64 






- 182.6 




[GK20011 30 




12.48 


13.33 


14.57 


18.11 


— 92.3 


2.1 


Sz 33 


8.22 


9.10 


9.84 


11.19 


13.99 


- 9.6 


2.2 


[GK2001] 31 




12.75 


13.73 


14.84 


17.12 


- 68.7 


1.8 


ISO-Chal 237 


7.69 


8.55 


9.34 


10.75 


13.90 


< - 3 




[GK2001] 53 




13.73 


14.57 






- 15.7 


6.2 



Kerrvon^^Gomez 



Gomcz^^Cenvan 



Cambresv et al. 



2001) 
2001) 



1998) 



plates and thus optical spectra can be obtained. The existence of an optical object at the 
position of the infrared discovered source was the selection criterion used to define the 
present sample. 

[SGR2003] 1 is the only object not selected as a candidate very low mass member of 
the cloud by previous surveys. This object was observed simultaneously with [GK2001] 
30 (ISO-Chal 220). The central slit of the MOS was placed on [GK2001] 30 whereas the 
rest of the slits on optical stars in the field. One of these objects (identified as [SGR2003] 
1) turned out to be a new low mass cloud member judging from the strong emission in 
Ha (equivalent width = — 18 2.6). Table 1 provides th e coordinates of this object and 
Table 2 K and H magnitudes jGomez fc Kenvon IboOljl . [SGR2003] 1 was not detected 
at J and thus not placed on the color-color (J — H vs. H — K) diagram to determine 
likely near-infrared color ex cesses. Similarly Sz 33, a previously known T Tauri star of 
the cloud (see, for example, lLawson et al. Ill996|) lying in the close vicinity of [GK2001] 



31 (ISO-Chal 225), was observed on the same MOS frame as this object. 



3.2. The Spectra 

We mention that for five objects in our sample (ISO-Chal 126, [GK2001] 30, Sz 
33, [GK20011 31. and ISO-Chal 237) infrared sp ectra are available in the literature 



1 Gomez fc Persi 



2002 



Oomez fc Mardones 



2003). However most of these spectra are 



the Anglo-Australian Observatory. The plates were processed into the present compressed 
digital form with their permission. The Digitized Sky Survey was produced at the Space 
Telescope Science Institute under US Government grant NAG W-2166. Copyright (c) 1993, 
1994, Association of Universities of Research in Astronomy, Inc. All right reserved. 
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almost featureless at the spectral resolution used, thus not allowing a full characteri- 
zation of the objects. [GK2001] 30 shows Br7 in emission while Sz 33 have Pa/3, Pa5 
and Br7 in emission. The spectral slope of these sources, with exception of Sz 33, in- 
creases with the wavelength suggesting the presence of cool material probably in a cir- 
cumstellar disk. Optical low resolution spectra of these objects provide additional useful 
information to support their pre-main sequence status. In particular Ha in emission is 
a distinctive feature of T Tauri stars and of previously known young BDs. Furthermore, 
spectral classification schemes based on optical features are in general better developed 
than near-infrared systems and thus can provide more reliable estimates for the low-mass 
candidates. 

Figs. 1-5 show the spectra of the observed objects. Photospheric features such as 
molecular bands of CaH, TiO, and VO, as well as atomic lines of Na I and K I are 
indicated. These molecul ar and atomic features are common among M-type type stars 



I Kir knatrick et al. 



1991). 



In the following sections we will use the detection of specific spectral features (such 
as Ha in emission and the Na I absorption doublet) to unveil the nature of our sample 
objects. The strength of the TiO molecular bands at 6180 and 7100 A and the PC3 
index of iMartin et al. I 1)19991) will be applied to derive spectral types for these objects. 



3.2.1. Spectral Features 



Table 2 lists Ha equivalent width measurements for the observed objects. We estimate 
an uncertainty of the order of 3-5% due to noise in our spectra and errors in placing the 
continuum level. Five objects in our sample (ISO-Chal 126, [KG2001] 102, [SGR2003] 1, 
[GK2001] 30 and [GK2001] 31) show conspicuous Ha emission line with equivalent widths 
> 50A. Sz 33 and [GK2001] 53 have moderate emission in this line. For other two objects 
in the sample, [GK2001]8 and ISO-Chal 237, we do not see a strong Ha in emission and 
are able to measure only an upper limit to its equivalent widths. Whereas higher signal to 
noise and resolution spectra are required to check in detail for the presence (absence) of 
Ha in emission, these two objects are most likely foreground dwarfs or background giants. 
As a group, spectra in Figs. 1-5 show similar characteristics as M-type T Tauri stars and 



previously know young BDs in star 



Briceno et al. 



1998: 



forming regions (see, for e xample 



Comeron et al. 



2000: 



Martin et al. 



Sch wartz 



1977 



200JJ). In parti cular, we note 



that the equivalent widths of Ha are comparable with those measured by 



Oomeron et al. 



I 2000) for their sample of very low mass stars and BDs in Chamaeleon I. 



Martin et al. 



1 1996) in their study of BDs in the Pleiades found that, due to lower 
gravity, the Na I doublet (8183-8195 A) tends to be weaker in young sub-stellar objects 
than in field stars of identical spectral type. These authors analyzed objects with spectral 



types in the M3.5 - M9.5 range. Table 2 lists the Na I doublet equivalent widths for the 
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objects in our sample for which the red part of the spectrum is available and of good 
quality. We estimate an uncertainty of ~ 0.5 A in our measurements. Using the spectral 
types derived in the next section (see Table 3), we compare the Na I doublet equ i valent 



Martin et al. 



(1996, 



widths of our objects with those of Pleiades and field stars listed by 
-their Table 5). 

Three of the objects ([KG2001] 102, [GK2001] 30 and [GK2001] 31) have very modest 
Na I doublet equivalent widths, smaller than those of Pleiades and field stars of the same 
spectral type, confirming that they are very low mass young objects. The Na I doublet 
equivalent width of [GK2001] 53 is comparable to those of two Pleiades BDs of the same 
spectral type, while the equivalent width of [GK2001] 8 is larger than that of the Pleiades 
BD Calar 4 (with identical spectral type as [GK2001] 8) and comparable to that of field 
dwarfs. This supports the hypothesis that [GK2001] 8 is not a low mass member of 
Chamaeleon I and is likely a f oreground dwa r f. Fin ally, Sz 33 with a M2 spectral type is 



outside the range analyzed by 



Martin et al. 



(1996). 



Finally, whereas the Li I A 6707.8 line is a well-known indica tor of youth in late- 
type fully convective stars (see, for example, 



Maerazzu et al. 



1993), the modest spectral 



resolution of our spectra, in combination with the low S/N ratio in a few cases, prevented 
us from detecting this line in our sample objects. Higher resolution data are thus required 
to unambiguously detect this line in the spectra of the present sample. 

In summary, Hcv together with the Na I doublet confirm that [KG2001]102, [GK20001] 
30 and [GK2001] 31 are very likely pre-main sequence, low mass cloud members, while the 
spectra of [GK2001] 8 and ISO-Chal 237 do not support cloud membership; the status 
of [GK2001] 53 remains uncertain, since it has an Ha comparable to that of previously 
known low mass cloud members, but a somewhat high Na I doublet equivalent width. 
However the equivalent width of this doublet for [GK2001] 53 is still comparable to those 
of two Pleiades BDs of the same spectral type. Finally, whereas we were not able to 
use the Na I doublet for Sz 33, ISO-Chal 126, and [SGR2003] 1, the presence of Ha 
in emission in their spectra supports the pre-main sequence natur e of these sources. I n 
addition Sz 33 is a previously known T Tauri member of the cloud ( Lawson et al. lllflfln ). 



3.2.2. Spectral Types 



In order to derive spectral types for our sources, we used the two TiO absor ption bands 



at 61 80 and 7100 A that are temperature se nsitive for M dwarfs and giants 1 0'Connell 
1973(1 . We also applied the PC3 index form iMartin et al. I <|1999(1 . For stars with avail- 



able flux-calibrated FORS spectr a, we compared th e m wi t h the detailed grid o f late 
type standard sta r s obt ained by 



Kirkpatrick et al. 



(1995). 



Kirkpatrick et al. 



Henrv et al. 



(1994), and 



8 
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Fig. 1. Optical spectra of [GK2001] 8 and ISO-Chal 126 obtained with the VLT+MOS 
and the CASLEO+REOSC, respectively. Molecular and atomic features common in M- 
type stars are indicated. The [GK2001] 8 spectrum has been flux-calibrated and normal- 
ized to the flux at 6800 A. 



Kenvon et al. 



1 1998) applied the Titanium Oxide band indexes at 6180 and 7100 
A to estimate spectral types for pre-main sequence stars in the Taurus- Auriga molecular 
cloud. Each index measures the depth of the TiO band at a given wavelength relative to 



an interpolated nearby continuum point. These authors defined [TiO]i and [TiO] 2 as 

F6180 



[TiO]i = -2.5 log 
[TiO] 2 = -2.5 log 



F6125 + 0.225(F6370 — F6125 
F7100 

F7025 + 0.2(F74qo — F7025) 



(1) 

(2) 
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Fig. 2. Optical spectra of [KS2001] 102 and [SGR2003] 1 obtained with the VLT+MOS. 
Molecular and atomic features common in M-type stars are indicated. The spectra are 
flux-calibrated and normalized to the flux at 6800 A. 



O'Connell lll973 K The 30 A bandpass filters used in these indexes are free from 



(see also 

strong emission lines (such as Ha) or telluric absorptions. In F ig. 6 we plot these indexes 



O'Connell 



(1973). Three 



for our sources along with the dwarf and giant sequences from 
objects (ISO-Chal 126, [SGR2003] 1 and [GK2001] 30) are not shown in the figure due 
to the large veiling effect that mainly affects the [TiO]i index as we discuss below. 

The new objects, on average, lie slightly above the sequences defined by the dwarfs and 



Kenvon et al. 



giants , having weak [TiO]i index relatively to the [TiO] 2 absorption band 
1 1998) noticed the same tendency in many T Tauri stars with strong emission lines 
in the Taurus cloud. The emission from a hot blue continuum source (with T ~ 10 4 ; 



10 



Saffe et al.: Cha I selected young stellar objects 



TiO CaH TiO VO KI 



35 
30 
25 

1 20 
^ 15 
xT 10 



[GK2001] 30 





6000 



7000 8000 
Wavelength in Angstroms 



9000 




-1 



TiO TIO CaH TiO VO KI 

n n i i 



Sz 33 





6000 



7000 8000 
Wavelength in Angstroms 



9000 



Fig. 3. Optical spectra of [GK2001] 30 and Sz 33 obtained with the VLT+MOS. 
Molecular and atomic features common in M-type stars are indicated. The spectra are 
flux-calibrated and normalized to the flux at 6800 A. 



Hartigan et al. 



1991) may diminish the depth of the TiO bands and induce us to assign 



earlier spectral types to the central star. This veiling effect increases toward shorter 
(bluer) wavelengths, affecting more the band at 6180 A than the index at 7100 A. As 
mentioned, five of the stars in Table 2 have Ha equivalent widths > 50 A which suggests 
the presence of veiling in their spectra. 

As to the luminosity class, pre-main sequence stars are known to have surface gravities 
inte rmediate between those o f dwarfs (luminosity class V) and giants (luminosity class 
III). 



Kirknatrick et al. 



1 1991) identified several features sensitive to the gravity of M type 



stars. In particular the intensity of the CaH absorption band at 6975 A and the Na I lines 
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Fig. 4. Optical spectra of [GK2001] 31 and ISO-Chal 237 obtained with the VLT+MOS 
and the CASLEO+REOSC, respectively. Molecular and atomic features common in M- 
type stars are indicated. The [GK2001] 31 spectrum has been flux-calibrated and nor- 
malized to the flux at 6800 A. 



at 8183-8195 A decrease with decreasing gravity and then are stronger in dwarfs than in 
giants. They defined two color ratios (A and C, respectively) that compare the fluxes in a 
30 A band centered on these features with respect to a nearby continuum integrated over 
an identical bandpass. These color ratios for our sample give values intermediate between 
those of dwarfs and giants and, on average, slightly closer to the giant branch. We caution 
that the veiling effect may also affect the depth of the gravity sensitive features and hence 
lead us to derive lower gravities, closer to class III luminosity objects. 
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TiO CaH TiO VO KI 

n r 



8 
6 
4 
2 

-2 



- 




[GK2001] 53 

u ill! 1 


1 




r 


6000 


7000 


8000 9000 



Wavelength in Angstroms 



Fig. 5. Optical spectrum of [GK2001] 53 obtained with the VLT+MOS. Molecular and 
atomic features common in M-type stars are indicated. The spectrum has been flux- 
calibrated and normalized to the flux at 6800 A. 



Three objects need to be discussed in detail. ISO-Chal 126 shows a featureless spec- 
trum (except for Ha and H/3, both in e mission, see Fig. 1 lower pan el); for this object we 



Gomez fc Mardones 



I 2003) from a near-infrared 



adopted a M2 spectral type derived by 
spectrum. 

The R-band spectrum of [SGR2003] 1 has a strong Ha line in emission and no other 
spectral feature and the TiO band at 6180 A is completely veiled (see Fig. 2 lower 
panel). The I-band region (not shown in Fig. 2 lower panel) was observed under poor sky 
conditions and some night sky lines were registered on the stellar spectrum. However, it 
clearly shows the TiO absorption at 7100 A, with a depth corresponding to a M6 spectral 
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13 




Fig. 6. [TiO]2 vs. [TiOli indexes for field dwarfs (filled circles) and giants (open circles) 



obtained from 



O'Cormell 



1 19731) . The continuous line connects the points defining the 
dwarf sequence and the dashed line those corresponding to the giant branch. Six of the 
young stellar objects in Tables 1 and 2 are indicated with starred symbols. 



type. Due to the veiling effect, it is likely that the depth of this band is diminished and 
hence, as discussed, we are assigning a spectral type earlier than the actual type. The 
I-band region also shows th e VO band centered at 7445 A that appears for spectral 
types as late as M7 or later i Kirkoatrick et ai~lll99fi ): for the latest M-sub types the VO 
band rather than the TiO band is indeed a more reliable indicator of temperature as the 
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1991). 



TiO starts to saturate IjKirkpatrick et al 
VO ratio as: 

0.5625(F 735 o - F 740 o) + 0.4375(F 7510 



Kirkpatrick et al. 



VO ratio 



7560 



1995) denned the 



(3) 



F7420 — F7470 

F\ 2 indicates the flux integrated between Ai and A2. For M7 stars they 



where F\ } 

obtain 1.02 < VO ratio < 1.07 and 1.07 < VO ratio < 1.12 for M8 stars. We derive 
VO ratio = 1.08 for [SGR2003] 1 indicating an M7-8 spectral type. We adopt an M7 
spectral type for this object; we caution however that we measured the VO band in a 
spectrum with sky emission lines superimposed on the stellar continuum. Although, in 
first approximation, this has a negligible effect on the depth of the detected absorption 
bands, better quality data are required to confirm our initial VO ratio measurement. We 
mention that a spectral type earlier than M7 would shift horizontally to the left this 
object on the HR diagram (see Fig. 7), corresponding to larger values of mass and age. 
In particular for an M6 spectral type we would obtain an age > 10 7 yr, which is hard to 
conciliate with the strong Ha emission and highly veiled spectrum of this object. Both 
features indicate that the object is most likely very young and probably not older than 
a few million years. In ad dition, the compari s on of the I-band spectrum of [SGR2003] 1 



Kirkpatrick et al. 



1 1991) indicates a spectral type M7 for this 



with the grid of spectra of 
source. 

In the case of [GK2001] 30 strong sky lines contaminated the region around the 
[TiO]i band and prevented us from measuring this index. In addition the R-band data 
are featureless (except for the presence of Ha in emission) indicating an increasing veiling 
effect towards this spectral region. From the depth of the TiO band at 7100 A we estimate 
a s pectral type M5 for t his object. 



Martin et al. 



1 1999) used the spectral index PC3 to derive spectral types for late- 
M and L field dwarfs. This index compares the fluxes in two 40 A bands centered at 
8250 and 7560 A, respectively. For spectral types M2.5 - LI they obtained the following 
relation: 



SpT = -6.685 + 11.715 x (PC3) - 2.024 x (PC3) 2 



(4) 



Martin et al. 



( 19990 call 



We calculated this index for our sample objects and used 
bration to estimate spectral types. 

Table 3 lists the spectral types base on the TiO-band indexes corresponding to both 
the dwarf (left side) and giant (right sid e) sequences, as well as those obtained from 



the index PC3 from 
the 



Martin et al. 



Kirkpatrick et al. 



1 1999) and those derived from the comparison with 
1 1991) 's grid. We also indicate the adopted spectral type for each 
object. For ISO-Chal 126 the near-infrared spectrum derived spectral type is given. In 
the case of [SGR2003] 1 we used the VO band rather than the TiO band to estimate 
spectral type as mentioned above. We estimate an uncertainty of about ± 1 subclass 



Saffe et al.: Cha I selected young stellar objects 15 

for the derived spectral types based on measurement errors of the TiO and VO bands. 
The same uncertainty is expected from spectral types deriv ed by the com p arison with 

(1999). The 



Kirkpatrick et al. 



Martin et al. 



1 1991) 's grid and from the PC3 index of 
veiling effect may increase this uncertainty up to ± 2 for the highly veiled and mostly 
featureless objects. 

As already mentioned, five of the objects (ISO-Chal 126, [GK2001] 30, Sz 33, 
[GK2001] 31, and ISO -Chal 237) have near-infrared spectra published a nd derived spec- 
tral types and masses IjGomez fc Mardones 20o3lG6mez fc Persi I200I . Excluding ISO- 



Chal 126 for which we adopted the near-infrared spectral type, for [GK2001] 30, ISO-Chal 
237, and [GK2001] 31, the optical and near-infrared derived spectral types agree within ± 
1 sub type. For Sz 33 we inste ad estimated a MO spectral type from our optical spectrum, 
while iGomez fc Persi I (|2002|) obtained a M2 type from the measured near-infrared water 
vapor bands. 

In order to derive effective temperatures, intrinsic colors, and bolometric corrections 
for our targets from the estimated spectral types, we adopted the calibrations obtained 



by 



Wil king et al. 



1 1999). These relations are valid for spectral types from M2 to M9. 



4. Luminosities, Masses and Ages 

We used the J magnitudes to estimate bolometric luminosities. This spectral band is less 
affected by contamination from circumstellar infrared excess emission than the K or H 
data. The bolometric luminosities are calculated from the following expressions: 



Log(L bol /L ) = 1.89-0.4M b oi 



M bol = J - Aj - DM + BCj 



Aj = 2.63E(J - H) 



(5) 
(6) 
(7) 



where DM = 6.0 is the distance modulus (d=160 pc: IWhittet et al. Ill997ft . BCj is the J- 
band bolo metric correction, and J i s the apparent magnitude. Aj is the extinction derived 
from the 



Rickc & Lcbofskv 



1 19851) reddening law, with E(J - H) = (J - H) - (J - H) Q 
the color excess, and (J — H) Q the intrinsic color. For two objects, [SGR2003] 1 and 
[GK2001] 53, wit h no J magnitude measur ements we used the K band data in com- 



bination with the 



Rick e fc Lebofskv 



(1985J) reddening law, A K = 1.79E(H-K) (with 
E(H - K) = (H - K) - (H - K) and (H - K) Q the intrinsic color) and the K-band bolo- 
metric correction, BCk, to estimate the bolometric luminosities. 

Luminosities for three of the o bjects (ISO-Cha I 126, [GK2001] 30, and [GK2001] 
31) were previously determined by 



Persi et al. 



12000). Our determinations agree with 



those of Persi et al. within a factor of ~ 2, with exception o f ISO-Cha I 126. F or this 
object, the main discrepancy comes from the estimated Aj. 



Persi et al. 



(2000) based 
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Table 3. Derived Parameters for the Observed Objects 



Name 


S.T(TiO) a 


S.T(PC3) b 


S.T(K91 grid) c 


S.T(adopted) 


Teff e 


A/ 


log Lboi/LQ 


M/M 


A 


ge 


[yr] 


[GK2001] 8 


M7/M4 


M5 


M4.5 


M5 
















ISO-Chal 126 








M2 d 


3426 


3.4 


0.34 


0.17 


< 1 


X 


10 s 


[KG2001] 102 


M8/M6 


M6 


M6.5 


M6 


2762 


1.1 


-1.071 


0.09 


4 


X 


10 5 


[SGR2003] 1 




M7 


M7 


M7 


2592 




-2.49 s 


0.019 


2 


X 


If)* 1 


[GK20011 30 




M6 


M5.5 


M5.5 


2845 


1.7 


—1.69 


0.075 


4 


x 


10 6 


Sz 33 


M2/M2 




M2 


M2 


3426 


1.8 


-0.12 


0.22 


1 


X 


10 s 


[GK2001] 31 


M2/M3 


M3.5 


M3 


M3 


3260 


1.4 


-1.81 


0.22 


4 


X 


10 7 


ISO-Char 237 


MO/MO 






MO 
















[GK2001] 53 


M6.5/M5 


M5.5 


M5 


M5.5 


2845 




-2.02 s 


0.06 


6 


X 


10 6 



a Left and right s ides from dwarf and giant branches, respectively, 
b Estimated from 



c Estimated from 



Martin et al. 



1999) 



Kirkoatrick et al. 
Gomez fc Mardones 1 1 2003^ 



Obtained from 



O'Connell 



(1973) 



( 199lh 's grid 



Wilking et al. 



(1999) calibration 



f We obtain A K = 1.0 and 0.895 for [SGR2003] 1 and [GK2001] 53, respec tively 
s Derived using K and H magnitudes and Ak = 1.79E(H — K) obtained from 



Rieke fc Lebofskv 



1985) 



their determination on the (I — J) color obtained by 
we use the (J — H) color inde x from 



4 and 3 in 



Persi et al. 



Cambresv et al. 



Gomez fe Kenvon 



(1998) while here 



(2001) to compute Aj (cf. Tables 



2000). Additional measurements of the I and H magnitudes may 



help to clarify this apparent lack of agreement. 

We chose pre-ma in sequence evolutionary tracks and isochrones from 

1997J) 4 to derive masses and ages for the observed sources. 



D'Antona fc Mazzitelli 



Gomez fc Mardones 



2003) used these models to study the m ass and age distributions 



of the young stars belonging to the Chamaeleon I dark cloud. 



D'Antona fc Mazzitelli 



1 1997) 's calculations provide consistent and plausible results for the group of objects 



they analyzed in reason able agreement with the higher mass members of the cloud 
(see 



Lawson et, al. 



1996). Table 3 lists these parameters. Fig. 7 shows the locations of 



observed Chamaeleon I stars on the HR diagram. 

We excluded from our analysis [GK2001] 8 and ISO-Chal 237 as they are probably 
not young cloud members based on the Ha emission and Na I intensity in the case of 
[GK2001] 8 (see Sect. 3.2.1.). Additional data and in particular higher resolution spectra 
around Ha are necessary to establish the true nature of these objects. On the other hand, 
we include [GK2001] 53 in this section although the Ha and Na I doublet equivalent width 



4 Available at http://www.mporzio.astro.it/~dantona/prems.html 
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3.6 3.5 

Log T eff 



Fig. 7. HR diagram showing the positions of the seven of the observed stars (starred sym- 
bols). Pre-main sequence evo lutionary tracks, indicated with continuous lines, are from 
D'Antona fc Mazzitelli I 1997L The thick continuous line corresponds to 0.08 M© (~ 80 



Mj), the H burning limit. The dashed lines correspond to the isochrones calculated also 
by the same authors. The horizontal error bars correspond to a ± 1 subtype uncertainty 
in the assigned spectral types. 



measurements do not allow us to unambiguously assert (disregard) the pre-main sequence 
nature of this object. 

Four objects in our sample ([KG2001] 102, [GK2001] 30, [SGR2003] 1, and [GK2001] 
53) have masses close or below the H-burning limit (i.e, 0.08 M©, ~ 80 Mj). In Fig. 7, the 
horizontal error bars indicate the displacements on the HR diagram that these objects 
would experiment due to ± 1 sub-type (i.e., AT ~ 150 K) uncertainty. [SGR2003] 1 is 
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a sub-stellar mass object regardless of the spectral type (or temperature) uncertainty. 
[KG2001] 102, [GK2001] 30, and [GK2001] 53 have masses close to the border line between 
the stellar and sub-stellar regimes. An uncertainty of ± 1 sub-type prevent us from 
establishing the (stellar or sub-stellar) nature of the objects (see Fig. 7). The rest of the 
sources have masses between 0.22 and 0.17 MQ. 

Luminosity determinations have a "typical" uncertainty of about a factor of 2 (or ~ 
0.03 in log scale). However the effect on our mass determinations is less significant than 
the uncertainty in spectral type or effective temperature as the evolutionary tracks are 
almost vertical. 

The masses for the objects for which near-infrared spectra are available agree within 
a factor of ~ 2 with those derived from the spectral classification based on near-infrared 
spectra, with the differences in the inferred masses reflecting the differences in the as- 
signed spectral types. 

Our sample sources span a range of ages betw een ~ 10 5 yr and fe w x 10 7 yr, in 



agreement with higher mass members of the cloud IjLawson et al. 



1996). In particular, 



only a small fracti on of our objects (two o ut o f nine) have ages > 10 7 year s, supporting 



the conclusions of 



Comeron et al. 



and 



Gomez fc Mardones 



on the star 



formation rate in the cloud. 

For the four objects in Table 3 with masses close or below the H-burning limit we 
compare our mass and age determinations, based on the 



D'Antona fc Mazzitelli 



(1997) 



models, with two addition al sets of pre-ma i n sequ ence e volutionary tracks and i sochrones 



We select the models of 



Burrows et al. 



Chabrier et al. 



(2002) 6 as both 



H997h 5 and 

include detailed calculations below the 0.08 Mq limit. Fig. 8 shows these tracks and 
isochrones superimposed on the HR diagram. In Table 4 we compare masses and ages 
derived using these models. [SGR2003] 1, [GK2001] 30, and [GK2001] 53 have sub-stellar 
masses, independently of the adopted calculations, while [KG2001] 102 may be stellar or 
sub-stellar depending on the tracks used. 



5. Summary and conclusions 



We present optical spectra of eight candidate very low mass members and a previously 
known T Tauri star of the Chamaeleon I dark cloud. All but two of these objects have 
Ha in emission with equivalent widths simil ar to other known v e ry low mass T Tauri 



stars and previously identified young B Ds La wson et al 



Neuhauser fc Oomeron 



1998; 



Martin et al. 



1996; 



Bricefio et al. 



1998; 



20011) . indicating their young age. The Na 



I doublet (8183-8195 A) is weaker than in stellar objects of the same spectral type for 
[KG2001] 102, [GK2001] 30, and [GK2001] 31, supporting the very low mass nature of 



5 Available at http://zenith.as.arizona.edu/~burrows/dat-html/data/ 

6 Available at ftp://ftp.ens-lyon.fr/pub/users/CRAL/ibaraffe/ 
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Fig. 8. Pre-mai n sequence evol u tionar y tracks calculated by 



left panel- and 



Chabrier et al. 



Burrows et al. 



-right panel-, indicated with continuous lines, 
and superimposed on the HR diagram. The thick continuous line corresponds to 0.08 
M© (~ 80 Mj), the H burning limit. The dashed lines show the isochrones calculated 
also by the same authors. The positions of four of the observed stars (see Table 4) are 
indicated with starred symbols. 



these objects iMartm et al. In the case of GK2001] 53 the Na I doublet strength is 

comparable to two known BDs in the Pleiades. Higher resolution spectra of these objects 
would allow the detection of the Li I A 6707.8 and thus an independent confirmation of 
their youth. 

The lack of significant Ha emission in ISO-Chal 237 and [GK2001] 8 casts doubts 
on the nature of these objects. The Na I doublet strength for [GK2001] 8 suggests that 
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Table 4. Comparison of Masses and Ages Corresponding to Different Models 



Name 


M/Mq 




e [yrs] 


M/M 


Age [yrs] 


M/M 


Age [yrs] 




B97 a 






C02 b 




DM97 C 




[KG2001] 102 


0.027 




< 10 5 






0.09 


4 x 10 5 


[SGR2003] 1 


0.016 


5 


x 10 6 


0.025 


7 x 10 6 


0.019 


2 x 10 6 


[GK2001] 30 


0.045 


7 


x 10 5 


0.04 


1 x 10 6 


0.075 


4 x 10 6 


[GK2001] 53 


0.035 


4 


x 10 6 


0.04 


5 x 10 6 


0.06 


6 x 10 6 



' Burrows et al. 


1997 




Chabrier et al. 


2002 


) 


'DAntona & Mazzitelli 



this object is probably a field dwarf. Higher resolution spectra are required to properly 
classify both ISO-Chal 237 and [GK2001] 8. 

In summary, out of eight candidates, two turned out not to be cloud members, the 
nature of other object ([GK2001] 53) would require additional confirmation by higher 
resolution spectroscopy, while Ha, together with the Na I doublet strength in some cases, 
support the pre-main sequence nature and low mass classification of the remaining five 
sources. 

We used the stren gth of the TiO or VO mole cular bands as well as the PC3 index from 



Martin et al. 



1 1999|) a nd the grid of spectra of Kirkpat rick et al. 



tral types and adopted 



D'Antona fc Mazzitelli 



(1997) 



1 1991), to derive spec- 



evolutionary tracks and isochroncs 



to estimate masses and ages . We detect one new BP and three transition ste llar /sub 



stellar objects of the cloud. 



Burrows et al. 



(1997) and 



Chabrier et al. 



(2002) models 



confirm these results. The rest of the observed objects have masses between 0.22 and 
0.17 MO. [SGR2003] 1, with a derived mass of 0.019 M© (~ 19 M, T ), lies very close to 
the deuterium burning limit (~ 13 Mj). This is the least massive object identified in the 
cloud so far. These sources span a range of ages be tween ~ 10 5 yr and few x 10 7 yr, in 



agreement with higher mass members of the cloud IjLawson et al. 



1996). 



The Chamaeleon I dark cloud has been searche d at different waveleng t hs to detect 



1999; 


Persi et al. 


2000; 


Gomez & Kenvon 


2001 





Comeron et al. 



1998; 



Oasa et al. 



(2000) reported the 



detection of four bona-fide BDs in the Chamaeleon I d ark cloud in addition to eight 
candidate transition st el lar /sub stellar objects (see also 



Neuhauser et al. 



2002). 



Gomez fc Persi 



Neuhauser fe Comeron 



Gomez fc Mardones 



1998: 



(2003) iden- 



I2002J) and 

tified ~ 30 very low mass objects, including several transition objects. These in- 
vestigations have significa ntly increased the low mass population of the cloud (see 



Gomez fc Mardones 



2003). However, these authors attempted a determination of the 



IMF of the cloud and noticed an apparent scarcity of very low mass objects in the whole 
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Chamaeleon I cloud based on the beha vior of the IMF in the central 300 arcmin 2 region, 
determined by 



Comeron et al. 



12000). This small portion of the cloud has been inten- 



sively surve yed at different wavelengths combining both photometric and spectroscopic 



information. 



Gomez fe Mardones 



I 2003) estimated that about 100 objects remain to be 
identified in the whole cloud. 

The results presented in this paper contribute to the census of the pre-main sequence 
population of the cloud. In particular, we identified one sub-stellar object, increasing 
the number of bona-fide BDs in the cloud from 4 to 5. The newly discovered BD as 
well as the three transition objects analyzed in this contribution will eventually help to 
determine the behavior of the IMF for the complete cloud in the sub-stellar regimen and 
a comparison with other nearby star-forming clouds. 

The newly detected objects increase the overall number of young BDs and transition 
objects known today. These objects will eventually allow us to carry out statistical studies 
of the physical properties of young BDs. To better understand similarities and differences 
among stars, brown dwarfs and planets, confrontations of the characteristics of the three 
groups are needed. 
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